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In the system during the first order phase transition new supercritical 
formations of a new phase appear with some fixed probability, but appear 
in stochastic manner. In a system with macroscopic sizes due to a giant 
value of the Avogadro number there appears rather big number of droplets. 
It allows to use the averaged characteristics to construct kinetics of a nu- 
cleation process. Kinetics based on averaged characteristics is described in 
p. In this paper the time evolution is constructed and it is possible to 
extract elementary intervals where thermodynamic parameters and the nu- 
cleation rate have small variations. The total number of droplets is so big 
that at every elementary interval there appears a great number of droplets 
AN. On the base of traditional thermodynamics one can state that the rela- 
tive fluctuation 6 AN /AN of droplets formed at elementary interval is small 
and has an order of (AA^)^/^. This remark completely solves a problem of 
justification of nucleation description based on averaged characteristics. But 
recently there appears a set of papers [0, where a stochastic effects (the 
effects of fluctuations of droplets formation) were described and investigated. 
Two approaches were formulated there. Authors didn't hesitate that these 
approaches , gave different results. So, one has to analyze approaches 
[|2|, [H] and decide what is the true result. 

An idea formulated in 0, allows to establish corrections to the total 
number of droplets A^ appeared in the system. It was supposed that these 
corrections are functions of A^. To demonstrate the error of this approach it 
is sufficient to take two identical systems then to calculate them separately 
and to add results or to calculate correction directly for the total system. 
The results will differ. 
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One has to determine a real volume to which one has to refer the number 
of droplets. It is simple to do with the help of results from p| . In that paper 
kinetics of nucleation for diffusion regime of droplets growth was constructed. 
It was shown that a solitary droplet perturbes vapor up to distances of an 
order y/ADt, where D is diffusion coefficient, t is a time from droplet forma- 
tion. One can take as t a time ti of nucleation period duration. A nucleation 
period is a period of intensive formation of droplets. 

It allows to give a new definition of the volume Vi where the number Ni 
of droplets is formed. Namely this value has to be regarded as a volume of a 
system. This volume is 47r(4i5i(:i)^/^/3. If the sizes of the system are smaller 
than this value one has to take the volume of the system as this value. But 
such a small system can be hardly regarded as a macroscopic one. At least 
one has to analyze carefully the boundary conditions. 

Naturally, the droplets appeared at different times perturbe initial phase 
at different distances. So, one can regard above formulas only as estimates. 
Some more rigorous equations can be found in |Q. 

The number of droplets A^^i isn't too big as is. So, an analysis of 
stochastic effects has a real sense. It is interesting now to get all correction 
terms which are ascending with the number of droplets (but not only a leading 
term). To solve this task one has to modify approaches from [|l|, p| where 
only a vague conclusion about the order of the leading term was made. 

Complexity of this problem lies in conclusion that one can not directly use 
equations based on the theory with averaged characteristics. In [Q, some 
properties of solution on the base of averaged characteristics were starting 
points for constructions. This supposition was adopted without any justifi- 
cations. 

We shall consider the situation of decay. The new dimensionless parame- 
ter - the number of droplets destroys the universality observed in ||^ for the 
theory based on averaged characteristics. Moreover, it is difficult even to 
formulate the system of equations. It radically complicates the problem. 

The property of effective monodispersity formulated in was used in 
without any justifications. Generally speaking this property can not be 
directly used to calculate stochastic corrections. This leads to an error made 
in [0. Now we shall formulate more correct approach. 

Both approaches from [0 and used the following property "The droplets 
formed at the beginning of the nucleation period are the main consumers of 
vapor" . This property is valid, but it is substituted in , [0 by the following 
statement: "The main source of stochastic effects are the free fluctuations of 
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droplets formed at the beginning of the nucleation period, they govern the 
fluctuations of all other droplets". But the last statement isn't valid. So, one 
has to use some new methods which are presented below. 

The use of monodisperse approximation will lead to some errors. But due 
to universality of solution these errors can not be be too big. Qualitatively 
everything is suitable, but precision will be very low. 

The same conclusion will be valid for calculations based on some model 
behavior of supersaturation (justification is valid for a vapor consumption, 
but not for stochastic effects). Here the final result will be more precise but 
it comes from rather spontaneous artificial choice of some parameter which 
equals in |l| to 1/2. In it is supposed that until some moment (it is 
chosen in [|T| as a half of nucleation period) the droplets are formed under 
ideal conditions and namely these droplets determine a vapor consumption. 
This approach taken from was used in in slightly another sense. It 
is supposed that droplets formed during a first half of nucleation period are 
the main source of stochastic effects. The last statement was not justified 
in and it is rather approximate. The relative correctness of a result was 
attained due to specific compensation of different errors of approximations 
used in 

All arguments listed above lead to necessity of reconsideration which will 
be made in this publication. A plan will be the following 

• On the basis of algebraic approach we shall see that stochastic effects 
are small 

• A smallness of stochastic effects allows to seek the solution on the 
base of the theory with averaged characteristics. But we have to take 
stochastic effects from all droplets formed during the nucleation period. 

• The possibility to take into account the role of all droplets can be 
ensured by the property of similarity of nucleation conditions during 
the nucleation period. This property can be considered in two senses - 
1) in the local differential sense and 2) in the integral sense in frames 
of the first iteration. The local property will be used in justification 
of the smallness of stochastic effects and the integral property will be 
used to calculate the stochastic corrections. 

All analytical results will be checked by computer simulation and a coin- 
cidence will be shown 
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All mentioned constructions will be valid for an arbitrary first order phase 
transition. 

The law of droplets growth will be a free molecular one, the linear size 
grows with velocity independent from its value. Consideration of other 
regimes can be attained in frames of the current approach by some trivial 
substitutions, but one has to take into account that the new regime requires 
new approaches to construct nucleation kinetics as it is shown in P]. So, we 
can not agree with he statement in that one of results is an account of 
stochastic effects in a diffusion regime of droplets growth. This effect has to 
be taken into account by application of methods presented in 

Dynamic conditions can be easily considered by direct generalization of 
methods presented here and we needn't to present it in details. 

1 Estimates for stochastic effects 

We consider the theory based on averaged characteristics. It is supposed to 
be known 0, that the supersaturation ( behavior can be determined after 
certain renormalizations by the following equation 



A good approximation for solution and for a distribution exp{ip) over linear 
sizes is fi = exp(— 2;'^/4). The form of fi is given by fig.l It is seen that at 
zo = 1.25 the nucleation period stops. 



For any moment t or z a function g has one and the same form. This is 
a similarity of nucleation conditions. We see that every time the droplets 
formed at the last third of a period from beginning of nucleation until a 
current moment will accumulate a negligible quantity of substance. The 
relative quantity of the substance there has an order of (~ 1/27) and is so 
small that even if there will be fluctuations the quantity will be small. 

From the form of /i it is seen that until z= = 0.7 zq all droplets will deplete 
vapor rather weak. It will be important for future analysis. 




This approximation has rather high precision 0. 
following law of substance accumulation 



It is based on the 
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Figure 1: A form of approximation for size spectrum 

The mentioned property of g allows to use a monodisperse approximation 
P] not only at the end of nucleation but in every moment of the nucleation 
period [^]. Let t{G) be the moment when there are G molecules in droplets 
(in appropriate units). An application of the monodisperse approximation 
0] leads to 

-Niz)z^ = a 

4 

Now it will be possible to repeat all constructions with G instead of 
1 (in renormalized units, before renormalization it would be l/F (see 
The sense of these transformations is rather evident. Let t{z) be the current 
moment of time {z is the coordinate of ta spectrum front, actually t is pro- 
portional z). We suppose that before az (a is some parameter) droplets are 
formed without mutual influence and one can write Poisson's distribution. 
Then a natural restriction on a appeared, namely a < 0.7 * zq/z. Then we 
suppose that the influence of other droplets on its own formation is negligi- 
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ble (this follows from 2 * 0.7 = 1.4 > 1.25 and from notation about the last 
third of nucleation period). Then it is possible to write Poisson's distribution 
for the second group of droplets, but with parameters depended on stochas- 
tic values - characteristics of the droplets distribution from the first group. 
Rigorously speaking one has to use the first four moments of the droplets 
distribution in accordance with 0, but for simplicity we shall use here only 
the zero momentum. As a compensation for this simplicity we has to use 
here a = 0.25. 

Then one has to come from Poisson's distributions to Gauss distributions 
and integrate them with account of connection between stochastic parameters 
of embryos formation from the first group and parameters of distribution from 
the second group. Unfortunately, in ||^ the final result was an expression for 
distribution, but not for the number of droplets which can be observed in 
experiment. Beside this the expression for the distribution was calculated 
only in a leading term which is certainly a Gauss distribution. Corrections 
haven't been established. 

Contrary to we shall take into account all correction terms which 
comes from transition from Poisson's distributions to Gauss distributions 
and corrections for nonlinear connection between group distributions. We 
shall take all terms which are growing when the total number of droplets 
grows. We shall get the following result for droplets distribution 

„ . 9a ,-,,2 , 9a D'^.. , 



'2nN{3a + l)' 3a + 1 2 

where 

D 
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N - some stochastic value of the total number of droplets, N - the mean 
value of droplets and y is the correction for spectrum. 

At =1/4 we get 

y = -^DJ8087DI - 10269)s + (-- + 

^ 74088 ^ ' ' ^ 9 1176 ' 
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is a small parameter of decomposition. To get all ascending corrections we 
must fulfill decomposition until s^. At arbitrary a we get the following ex- 
pression 



1 + WiS + W2s'^ 



Here 



wi = -^((486Z^ + mDll^'^ - %m^^Dl + 324/^ - 648D^Z^° + 756/^^ 

+810Z^D^ - 459^^ + 2mDll^ - I2,hf D]- 
387f+m^Dl-153l^+27l^Dl-2m^-66l^Dl-3l^+27l^-3l^Dl+m^Dl-3l^+ 
91 + IDl + Dl)Ds)/{{l + 3/2)3(/ + 1)(-1 + 



where 
and 



I — y/a 

W2 = Wo2/{{l + 1)'(-1 + + 31^1') 

3 
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Having integrated this expression we get corrections to droplets number. The 

term at s gives zero after integration and the first correction has an order of 

~ and doesn't depend on the total number of droplets. A coefficient at 

has at a = 1/4 a value 

do = 311/3024 

At arbitrary a a coefficient at in correction for the total number of droplets 
will be 

4 = ;^[108a^ + 540a^^/2 - 72a^ - 9300*^/2 - 336a^+ 

713a^/2 + 158a^ + Aa^ - 2030^/^ _ ga + 390^/=^ - 3a^/^]/[a^/^{l - af{l + 3a)] 

It will be interesting to compare results with and without corrections from 
transition from Poisson's distribution to Gauss distribution. At the leading 
term there will be no change. At correction terms we have 

732037 ^4 24137569 . 13, 2 



12446784 ' 10978063488 ' 36' 
do = -37/126 

4 = -l(648a^^/2 - 216a^ - 10620^^/^ + 108a^ + 753 J/^- 

1 ^ 
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Figure 2: Dependence of corrections on a 

30a^ - 195a^/^ - 12a^ + 19a^/^ + 6a - 7a^/^)/{a^/^{l + 3a) (1 - a)^) 

It is seen that these corrections are small. At arbitrary a except too small 
and 1 (these values are unreal) we get values shown at fig. 2 

Points show corrections with transition from Poisson's distribution to 
Gauss distribution, a line shows corrections without transition from Poisson's 
distribution to Gauss distribution. Both corrections have one and the same 
order and they are small. 

One has to note that we have not take into account corrections from tran- 
sition from summation to integration. It is made due to reasons formulated 
below. Really, we have at transition from summation in formula 

P{N) = ^ Pi(7Vi, N,)P2{N - N,, N2{N,)) 

Ni 

to integration 

P{N) = J dAP,{N,,N,)P2{N - A,N2{N,)) 
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to use Euler's decomposition. But discrete character in nucleation isn't so 
trivial. The process of vapor consumption can not begin without the first 
droplet. The system will wait for droplet as long as it will be necessary. It 
shows that discrete effects are complicate and require a separate publication. 

To use Poisson's distribution for the first group one has to make the 
following notation. Really nucleation conditions for the first group don't 
differ from the whole group. So, for distribution for the first group one has 
to take distribution Pi with reduced halfwidth. But one can not attribute 
a halfwidth to Poisson's distribution. That's why we considered effects with 
and without corrections from transition from Poisson's to Gauss distribution. 
So, we can use Gauss distributions as initial ones. For Gauss distribution 
one can easily reconsider the halfwidth. Then for Pi one can take 

Pi ~ exp((7Vi - Ni)y{2bNi)) 

where 6 is a renormalization coefficient. Distribution P2 remains previous 

P2 ~ exp((7V2 - A^2)7(2A^2)) 

where N2 is given by 

Ar 1 2 2 14 3 35 4 91 5 

N2 ^ (1 - -s + -s^ H - a)N 

^ 3 9 81 243 729 ^ 

where N is a. mean total number of droplets, 

Ni-aN 



s — 



aN 



is a small parameter of an order N 
After integration one comes to 

P ~ expf —d ) 

2(9a + 6 - 66a + 9ba^ - 9a^) ' 

where 

d= {N -N)/^fN 

The halfwidth of the distribution Pi must be equal to the halfwidth of 
P, which leads to 

«(!-«) 



6 = 9- 



-9a2 + 15a - 1 
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Now we shall consider effects of renormalization. The ratio of corrections 
with renormalization and without renormalizations is given by 




Figure 3: Ratio of half widths 

For all reasonable values of a it is approximately 1. At a = 1/4 we get 
1.0041 So, here the effect of similarity of nucleation conditions doesn't lead 
to remarkable effects. It is only due to monodisperse approximation. Later 
this effect will be essential. 

Instead of taking into account all moments of distribution we can directly 
calculate the effects on the base of exphcit form of spectrum in frames of 
iteration procedure. 

The result of this section is the following: we have proved that stochastic 
effects are small. Beside this we have demonstrated how to use the similarity 
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of nucleation conditions. Here it is useless, but later this method will lead 
to some essential numerical results. 



2 Stochastic effects in the iteration proce- 
dure 

There is only one parameter which is essentially deviated by stochastic effects 
- it is dispersion of distribution or a halfwidth of distribution. Now we shall 
calculate this value. 

The most advanced approach was suggested in |^ . But even this approach 
has many disadvantages. 

We shall characterize a droplet by linear size p which is the cubic root of 
its molecules number. Its velocity of growth does not depend on p. 

Decomposition of a whole interval of nucleation into elementary intervals 
is connected with some difficulties. An elementary length A must satisfy 
according to IQ to two requirements: 1. A number of droplets formed during 
elementary length must be big. 2. An amplitude of a spectrum has to be 
approximately constant during an elementary interval. 

It is clear that the second requirement can not be satisfied. Stochastic 
deviations of an amplitude leads to violation of the second requirement. 

We shall apply the second requirement not to stochastic amplitude as it 
was stated in but to averaged amplitude. Then the second requirement 
is : An averaged amplitude of a spectrum has to be approximately constant 
during an elementary interval. 

Stochastic amplitudes fi are introduced in as 



where Ni is the number of droplets formed during A. It isn't stochastic 
value but partially averaged one. An expression for the number of molecules 
in droplets formed during interval number z at a moment tk (it means that 
now we are at interval number k) will be the following one 



The difference between forth powers corresponds to a constant amplitude of 
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spectrum. It is wrong and then eq. (12) in |]T|] and all further equations are 
not correct. 

But the is no necessity to use such a way to account the number of 
molecules in a new phase. It is absolutely sufficient to take the following 
expression 

which is valid at A; — i ^ 1. In a whole quantity of substance it is sufficient 
to take into account only droplets with k — i ^ 1. The relative weight of 
dismissed terms will be small. 

Then for the total number of molecules in droplets at interval number k 
we have the following expression 

i=l 

where Xk-i is appropriate coordinate or 

k 
i=l 

This representation is important because now the note in |^ after eq. 
(15) isn't necessary. This note stated that the probability for A^j to deviate 
from N is very low. This note is doubtful because namely these deviations 
are the base for stochastic effects. Now there is no need in this note. 

The next step is to build [0 two cycle construction for nucleation period. 
During the first cycle the main consumers of vapor appeared and during the 
second cycle they rule a process of formation of all other droplets. In it is 
supposed that during the first cycle a vapor depletion is negligible and during 
the second cycle new droplets are absolutely governed by droplets from the 
first cycle. Now we shall analyze an effectiveness of such procedure. 

We analyze an equation 

-0 = / (z — x)^ exp(—'ip))dx 
Jo 

The first iteration ^ is practically precise solution and it gives the number 
of droplets 

1 4y^7iV2 
Ntot = --r--V = 1-2818 
4 r 3/4 
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A model solution requires that until z = p there will be no depletion of 
vapor and then only droplets formed before z = p will consume vapor. Then 
for a total number of droplets we have an expression 




0.4 0.6 0.8 1 1.2 1.4 

P 



Figure 4: The ratio of mean numbers of droplets 

Now it is clear that in |l| the value of parameter of separation into two 
cycles was not chosen in a good style (at least from the point of view in the 
theory with averaged numbers). We can also stress the smooth dependence 

Ntot appr on p. 

We shall study the probability Pk of formation of stochastic number Nk 
of droplets at the first k elementary intervals 

Our constructions now resemble |1| but there is one essential difference. 
When we get a formula analogous to (30) we have no necessity to linearize 
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expression with respect to {Ni — Ni)/Ni, where Ni is a stochastic number 
of droplets formed at interval i, Ni is a mean number of droplets formed 
at interval i (it is a function of stochastic numbers of droplets at preceding 
intervals). This linearization can not take place because a ratio {Ni — Ni)/Ni 
can be zero or can attain huge value (with low probability). It is more simple 
and more justified to linearize expression on J2iPK^i " where pi is 

a linear size of droplets formed at interval i (all of them have approximately 
the same size). Really, due to summation variations of J^iPK^i ~ 
are much smaller than variations of [Ni — Ni)/Ni. 

Variations of {Ni — Ni)/Ni would be small only at very big numbers of 
droplets Ntot- One can get 

{N, - Ni)/N, ~ iV//' 
Ntot ~ MNi 

M is a number of elementary intervals. So, the theory with linearization 
proposed in [|l| would be well justified only in a region where the result can 
be obtained on the base of averaged characteristics. But the linearization 
proposed here leads to the same numerical expressions. For dispersion of the 
total distribution one can get 

DOC ^ ^~oon - ^) 

a 

where 

/3 = /3i-/32 
A = 8/ dU rfr(e-r)3exp(-e^) 

Jl/2 Jo 
/•oo POO /■1/2 

/?2 = 16 / dv dT{T - 0'(r - vfexpi-^^) exp( V) 

Jl/2 Jl/2 Jo 

/•oo 

a = (ixexp(— x^) 
Jo 

In two cycles construction the value of a, which is proportional to the 
total number of droplets has to be reconsidered and recalculated on the base 
of two cycles. Then we have to use instead of previous a = ao a new value 

/•oo 

a = ai = l/2+ / exp(-x^ + (x - 1/2)^) 

Jl/2 
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We shall use parameter p of separation of two cycles and we shall calculate 

«! as 

POO 

ai=p+ / exp(— + (x — pY) 
Jp 

Then according to fig. 4 we see that the ratio ao/ai is greater than 1 and ai 
is greater than ao- Here we see that two-cycles construction is approximate 
one. Result for D°° will differ from 0] and will be 

= iV°°0.69 

instead of 

^oo ^ JYOO0.67 

as it is stated in [|l|. The value differ from a real result more than one 
tenth. So, the new theory is necessary. 

Now it is necessary to decide what p shall we choose. At arbitrary p the 
expression for /3 will be the same but for /?i /?2 we have 

roc rp 

Jp Jo 

roo roo pp 

p2 = l6 di dri dT{T - if{T - 7]f exp(-e^) exp(-r/^) 
Jp Jp Jo 

We have to reconsider expression for a. 

After calculations we have for dispersion as function of k = p the following 

fig- 5 

A minimal dispersion will be at p = 0.6. This is the true value of p 
because this value corresponds to the sense of minimal work. Dispersion will 
be 

^oo ^ 2A^~°°0.66229 

We shall make our result more precise. We see that (5i and j32 are the first 
two terms of some series. We don't know other terms, but it is reasonable to 
suppose that this series resembles geometric progression with denominator 
P2/Pi- This leads to dispersion 

j^oo ^ 2A^~°°0.64107 

As it follows from fig. 6 the value of extremum remains p = O.Q 

An absence of extremum shift is important and is necessary for this ap- 
proach to be a self consistent. 



16 



0.26 



0.24 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 

k 



Figure 5: Dispersion as function of k 

Now we shall make the value of dispersion more accurate Due to the 
similarity of nucleation the first cycle doesn't differ from the whole period. 
Function j3 for the first cycle will be 

CP ppi 

A-S / dU rfr(e-r)3exp(-e') 
Jpi Jo 

rp rp rpi 

/32 = 16 / dU dv dr(r - 0'(t - exp(-e') exp(-77^) 

J pi J pi Jo 

Calculations for pi = 0.6 * 0.6 = 0.36 and p = 0.6 give p = 0.0255. 

Now we shall reconstruct (3eff based on Df" = 27V°°0.64107 We have 
(3^ff = (1 - D°°/N^)/a where a = 0.9092. Calculations give peff = 0.3268. 
One has to add f3 = 0.0255 which leads to (3 — 0.35156. For this value of /3 
the value of dispersion will be 

= 2Aroo0.61332 
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Figure 6: Dispersion as function of cl = p 



This is our final result. 

When we consider the dispersion which differs from — 2N°° we have 
the shift of result 

SD = Df - 

This shift has to lead to a variation of dispersion in a fist cycle and to a 
variation of the total dispersion. Having linearized this effect we get to the 
final value of dispersion 

r)oo 

4 2 ^oo 

and after calculations 

= 2A^~°°0.62309 
This value practically coincides with the previous value. 
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Another way to make results more precise is take into account the shift 
of dispersion directly in initial formulas. Having written for the dispersion 
correction in the first cycle 

with parameter 7, we get for the final distribution 



exp[- 



.[7VW-ArW-EjLi%(A^,-iV,)P 



2(A^W - PNi) 

where P is the number of elementary intervals until argument p, iV^*^) is the 
number of droplets calculated on the base of the theory with averaged char- 
acteristics, Ni the mean number of droplets formed during interval number i 
with account of fluctuations from previous intervals. The values are given 

j=p+i 

and M is the total number of intervals. 

Having fulfilled integration dNidN2...dNp, we get for a limit value of 
dispersion 

a a 

which leads to equation on 7, which can be easily solved 

1 - ^ 

7 



a a 



Calculations lead to 



7(p = 0.6) = 0.5083 

This result is very strange. It radically differs from the previous one. The 
reason of an error is that the duration of the first cycle is limited by p. So, 
we have to limit the duration of a whole period. The limit is, evidently, ~ 1. 
So, we have to recalculate /3 as 



19 



Jp Jo 

P2 = 16 /' [' dv r dr{r - 0'{r - r^f exp(-e^) exp(-77^) 
Jp Jp Jo 

Now we shall show a dependence Piniuai on p at fig. 7. 
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Figure 7: ^initial as a function of p 

Calculations give Pi„itiai — 0.17 and for the final dispersion 

= 2A^~°°0.6436 

This value coincides with a previous approach. It can be made more 
accurate by a geometric progression summation spoken above. 

3 Numerical results 

Numerical simulation of nucleation can be done by the following method. We 
split the nucleation interval into many parts (steps) . At every step a droplet 
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will be formed or not. The probability to appear must be rather low, then we 
have the smallness of probability to have two droplets at the same interval. 
This means that the interval is "elementary". 

The process of formation is simulated by a random generator in a range 
[0, 1]. If a generated number is smaller than a threshold parameter u, then 
there will be no formation of a droplet. If it is greater than a threshold, we 
shall form a droplet. As a result we have spectrum / of droplets sizes. Now 
it is a chain of and 1. The parameter u descends according to macroscopic 
law 

u = ■Uoexp(-rG/$) 

from a theory with averaged characteristics (it is based only on a conservation 
law without any averaging and can be used). Here 



$ is the initial supersaturation, Fc is a free energy of critical embryos forma- 
tion, i>c is a number of molecules inside a critical embryo, G is the number 
of molecules in a new phase taken in units of a molecules number density 
in a saturated vapor. By renormalization one can take away all parameters 
except G. 

To simplify calculations radically one can use the following representation 
g for G: 

G = z'^Gq — Sz'^Gi -\- 3zG2 — G3 
where z is a coordinate of a front of spectrum, and Gi are given by 

CZ 



Gi = fix)x^dx 
Jo 



/o 

We needn't to recalculate Gi, but can only ascend the region of integration, 
having added to integrals z^f{z)dx at every step. 

Our results are given below. The interval is split into 30000 parts. Pa- 
rameter Mo have been varied from up to 1 which leads to a different number 
of droplets. It is clear that the limit values are not good: at there are no 
droplets in the system, at 1 our intervals are not elementary. At every uo 
results were averaged over 1000 attempts. 

Shifts of droplets numbers are drawn at fig. 8 as a function of In A^(°°) 
It is seen that an analytical result about negligible value of corrections is 
correct. 
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Figure 8: Shifts of droplets numbers as a function of lnA^(°^) 

Dispersion as a function of In A^(°°) is shown at fig. 9. 

It is seen that analytical value of dispersion coincides with numerical 
simulation. The ends of the curve correspond to a zero number of droplets 
and to a giant number of droplets when the elementary intervals are not 
elementary and have to be thrown out. 

Stochastic effects in dynamic conditions @] can be analyzed by the same 
method. We needn't to describe it here. Numerical results are drawn below. 
Fig. 10 shows the shift in the number of droplets. It is small. Dispersion is 
drawn at fig. 11 (i.e. the value of 7). It is greater than in the case of decay. 

The physical reasons for the smallness of the droplet number shift for 
decay and for dynamic conditions will be different. 

For decay the reason is the following. The system wait the first droplet 
as long as necessary. Actually the time for kinetics of this system is t{G) 
with no connection with real time (certainly, the rate of nucleation has such 
connection). This phenomena is the reason for a smallness. 
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Figure 9: 7 as a function of In A^(°°) 

In dynamic conditions there is a time dependent parameter - the intensity 
of external source. So, there is no such a reason. 

But here in the theory with averaged characteristics there is a property 
of a weak dependence on microscopic corrections for a free energy 0. The 
same is vahd also for fluctuation deviations. So there will be a weak effect 
of stochastic nucleation. 

Because the reasons for smallness of effect in decay and dynamic condi- 
tions are different it is interesting to see whether they continue to act when 
supersaturation is stabilized at some moment. Analytical results shows that 
the will be an overlapping of two reasons. 

Really, if stabilization takes place at the period where the main consumers 
of vapor are going to appear then the majority of droplets appear in the 
situation when there is no influence on the system. Then the situation for 
these droplets resembles decay conditions (and may be even better because 
the external supersaturation is going to decrease). So the reason for the 
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Figure 10: Shift in a droplets number as a function of In A^(°°) for dynamic 
conditions 

decay situation works here. 

If stabihzation takes place at the second cycle, then the behavior of su- 
persaturation is governed by droplets formed in dynamic conditions and we 
have here the reason for smallness in dynamic conditions. In both situations 
the effect is small. Numerical results confirm this conclusion. 

The main result of this publication is a correct definition of all main 
characteristics of stochastic nucleation. It is shown that the main role in 
stochastic effects belongs to all droplets, but not to the main consumers of 
vapor. Only the property of the nucleation conditions similarity allows us to 
solve the problem of account of all influences during the nucleation period. 

When all disadvantages of , [0] are shown it is clear that these theories 
can not be considered as a solid base for nucleation investigation. 

But why results obtained in are so close numerically to real values. 
The reason is that on a level of averaged characteristics there is a universality 
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Figure 11: 7 as a function of In A^(°°) for dynamic conditions 

of nucleation process. So, the errors of [|1|, cannot lead to a qualitatively 
wrong results. 

One has to stress that all corrections obtained in this paper are also 
universal ones. Certainly, they are some coefficients in decompositions and 
the functional form is prescribed now (contrary to the theory with averaged 
characteristics). 

It seems that all effects considered here are negligible. For simple systems 
it is really true. But for systems with more complex kinetic behavior these 
effects can be giant. One of such systems is already described theoretically 
and this description will be presented soon in a separate publication. 

In diffusion regime of droplets growth one has to use another approach 
based on ||^, [0. In [|^ an explicit description of nucleation with account of 
stochastic effects was constructed. 

A nucleation with growing volumes of interaction will be presented in the 
next paper. 
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